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Using the exact Bethe ansatz solution of the Hubbard model and Luttinger liquid theory, we 
investigate the density profiles and collective modes of one-dimensional ultra-cold fermions confined 
in an optical lattice with a harmonic trapping potential. We determine a generic phase diagram 
in terms of a characteristic filling factor and a dimensionless coupling constant. The collective 
oscillations of the atomic mass density, a technique that is commonly used in experiments, provide 
a signature of the quantum phase transition from the metallic phase to the Mott-insulator phase. 
A detailed experimental implementation is proposed. 
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Introduction. — The Mott metal-insulator transition 
(MMIT) is a fundamental concept in strongly correlated 
many-body systems. Recent experiments with ultra- 
cold atomic gases in optical lattices are paving the way 
to explore such quantum phase transitions in a well- 
controlled manner. The MMIT with bosonic atoms has 
been demonstrated by Greiner et al A demonstra- 
tion with fermions has not yet been realized experimen- 
tally, although its realization is within reach of present- 
day techniques. In the fermionic case, the relevant theory 
is the widely studied Hubbard model. This is the sim- 
plest lattice model of interacting fermions, and is exactly 
soluble in one dimension. 

A gas of single component 40 K atoms in a parallel plane 
lattice has already been created by the LENS group, 
thanks to rapid progress in the cooling of fermions to 
temperatures below a micro-Kelvin. This demonstrated a 
non-interacting band insulator behavior 0, 0, Oil through 
the observation of dipole mode oscillations. An inter- 
acting gas of ultra-cold fermions with two populated hy- 
perfine levels in a one-dimensional lattice is also feasible, 
thus offering the possibility of observing the fermionic 
MMIT. 

Motivated by this opportunity, we address the prob- 
lem of how to detect the emergence of fermionic Mott- 
insulator phases in real experiments with ultra-cold 
fermions. With optical lattices, a harmonic potential is 
necessary to prevent the atoms from escaping, so that the 
Mott-insulator phase is restricted to an insulator domain 
at the center of the trap, and coexists with two compress- 
ible metallic wings. Therefore, the insulator phase can- 
not be characterized by a global compressibility as in the 
unconfined case. Rigol et al Q showed that a properly 
defined local compressibility exhibits critical behavior on 
approaching the metal-insulator boundary. 

In this Letter we show that collective oscillations of 
the atomic mass density, an indicator of compressibility, 
can be utilized to monitor the emergence of the Mott- 



insulator phase. We consider a zero temperature, one- 
dimensional Hubbard model with a harmonic potential, 
as a model of an ultra-cold two-component fermionic 
atomic cloud in a deep optical lattice with strong ra- 
dial and weak axial confinement. Based on the exact 
Bethe ansatz solution of the homogeneous ID Hubbard 
model |6], together with the local density approximation 
(LDA), we calculate the density profile of the cloud as 
functions of a characteristic filling factor and coupling 
constant. This leads to a generic phase diagram includ- 
ing a metallic phase and a Mott-insulator phase. We 
then investigate the low-energy collective density oscil- 
lations of the cloud in different phases using Luttinger 
liquid (LL) theory which describes long wavelength 
hydrodynamic behaviour. 

We find that in the metallic phase the collective oscilla- 
tion is an overall motion that goes through all sites of the 
cloud. This quenches gradually towards the phase transi- 
tion point, with the mode frequency decreasing monoton- 
ically to zero. After entering the Mott-insulator phase, 
the density oscillation revives, but is restricted to the 
compressible wings. Therefore, a sharp dip appears in 
all collective mode frequencies in the vicinity of the phase 
boundary, giving a clear signature of the MMIT. 

Phase diagram. — We consider the ID Hubbard model 
H, as a prototype for N (N/2 = = N±) interact- 
ing fermions in deep optical lattices, with a Hamiltonian 
given by: 



Ho/th = nil nil 
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Here th and u (u > 0) are the hopping parameter and the 
dimensionless on-site repulsion, respectively. The spin 
a =|, I represents different hyperfine states. The ax- 
ial harmonic trap adds a site dependent potential to the 
Hamiltonian so that: 
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where coo is the bare trapping frequency and d = A/2 is 
the lattice periodicity, while A is the wavelength of stand- 
ing waves used to create the ID optical lattice. We de- 
scribe the inhomogeneous gas, Eq. @, within the LDA. 
This amounts to determining the chemical potential of 
the gas from the local equilibrium condition, 



A* = /^hom [n 0) , u] + 



moo, 



°d 2 x 2 . 
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We also impose a normalization condition, 

as the chemical potential of a homogeneous system with 
a local filling factor n(x) (0 < n (x) < 2) [8]. In the 
framework of LDA, we replaced the site index i with a 
dimensionless variable x. The LDA is applicable in the 
limit of a large lattice with a slowly varying trapping 
potential. In other words, the use of LDA is justified 
if the size of the gas is much larger than the harmonic 
oscillator length ^ft/mujo, implying \i ^> c^o or TV ^> 1. 
Because of the optical lattice, the single-particle trap- 
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ping frequency is renormalized to u e ff = (m/m*) coq, 
with an effective mass m* = ft 2 / (2d 2 th)- We also 
define a dimensionless frequency, u = frw e ff/ (2^), with 
corresponding dimensionless time r = tu e ff/u. 

To construct the phase diagram, it is useful to intro- 
duce two parameters: 

• k = u 2 / (WNco) is the interaction strength 

• v = 2^/Nuu/tt is the characteristic filling factor 
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Figure 1: (color online) Phase diagram of one-dimensional 
confined Fermi gases in optical lattices. A: a pure metallic 
phase. B: a single Mott insulator domain at the center, ac- 
companied by two metallic wings. C: a metallic phase at the 
center surrounded with Mott insulator plateaus. D: a band 
insulator at the center with metallic wings outside. E: a band 
insulator at the trap center surrounded by metallic regions, 
in turn surrounded by Mott insulators. We have plotted a 
dashed line to illustrate the crossover behavior from A to C, 
or D to E. The empty symbols are the results from quantum 
Monte-Carlo simulations [5]. In Figs.(b), (c) and (d), the cen- 
tral densities n(x = 0) as a function of v have been plotted 
for three values of k — 0.01, 1, and 100. The dashed lines in 
these figures divide the system into different phases and the 
corresponding typical density profiles are shown in insets. 



We clarify the physical meanings of these parameters by 
considering two limiting cases. 

(i) In the limit of a dilute Fermi cloud [n(i) — > 
0], the Hubbard model reduces to Yang's exactly 
solvable Fermi gas model. Therefore, the homoge- 
neous chemical potential takes the form /ihom [n, u) = 
(2th) (7r 2 n 2 /8) /i h0 m [u/ (2n)] up to an irrelevant con- 
stant. The function ^' hom is obtained by solving a set 
of integral equations 01 . 

In terms of the dimensionless chemical potential fl = 
8/jl/ (u 2 th) and the rescaled coordinate x = Aoux/u, the 
local equilibrium and normalization conditions can be re- 
expressed as: 7r 2 / [27 2 (x)] l~i' horn [7 {x)\ + x 2 /2 = fit, and 
f+™2/j(x)dx = 16Nco/u 2 , where 7 (x) = u/[2n(x)}. 
These equations indicate that the coupling strength is 
k = u 2 1 (16iVa;), where k«1 corresponds to weak cou- 
plings and k ^> 1 is the strongly interacting regime. 

(ii) In the non-interacting limit, /ihom [n, u = 0] = 

-2t h cos (7ra/2) n S° 2t h (-1 + 7r 2 n 2 /8). It is straight- 
forward to sho w th at the central number density is 
n(x = 0) = 2\jN(jjji\, which we define as the charac- 
teristic filling factor v. 

We determine the phase diagram in Fig. la from the 
density profiles. Five phases (or mixed phases) are iden- 
tified by plotting the central density n(x = 0) as a func- 
tion of v for different values of as shown in Figs, lb, 



lc and Id, and explained in the figure caption. In addi- 
tion, various interaction regimes are illustrated explicitly 
by two dotted lines in Fig. la, where ft < 0.1 corre- 
sponds to weak couplings, n > 10, to strong couplings, 
and an intermediate regime is in between. In the weak 
coupling limit (ft = 0.01), the cloud is noninteracting. 
An increase of the characteristic filling factor changes it 
from a metal to a band insulator. In the strong coupling 
regime (ft = 100), the cloud behaves like a noninteract- 
ing spinless Tonks gas. It becomes a Mott insulator with 
increasing v, which can also be understood as a band 
insulator in a spinless gas. In the intermediate coupling 
regime, the existence of the trap leads to a coexistence 
of metallic and insulating phases, labelled C and E, that 
have no correspondence in the unconfined case. An ex- 
ample is shown in Fig. lc to illustrate the successive 
emergence of mixed phases C and E as v increases. 

A similar phase diagram of ultracold fermions in ID 
optical lattices has been proposed previously in Ref. 
using quantum Monte-Carlo calculations. Our results for 
the phase boundaries are in quantitative agreement with 
their findings in the same parameter space that has been 
simulated. This is expected since the LDA should work 
well for these parameters. Since our theory uses an exact 
analytic solution, together with the asymptotically exact 
LDA, this verifies the previous numerical work. The dis- 
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crepancy in determining the crossover boundary from A 
to C or D to E is due to the use of different criteria. 

The phase diagram in Ref. [5] covers only the weak 
coupling and intermediate regimes of n < 1.0. We have 
extended their phase diagram to the strongly interacting 
regime. This extension turns out to be crucial from the 
experimental point of view. As we shall see in Figs. 2a 
and 2c, the collective modes behave distinctly in differ- 
ent interaction regimes. Let us now consider in detail 
these collective density oscillations which are of primary 
interest in the present Letter. 

Collective density oscillations. — We describe the low- 
energy dynamics of density fluctuations of metallic re- 
gions in each phase using the LL model of the ID Fermi 
gas in optical lattices [7], where the Hamiltonian can be 
written as: 



dx 



u v {x) 



v=p,a 



1 



K v (x) \ dx 



(4) 

Here u p and u a are the density and spin velocities, re- 
spectively, and K p and K a are the Luttinger exponents of 
long- wavelength correlation functions. These parameters 
in a spatially homogeneous gas can be calculated from the 
Bethe ansatz solution of the ID Hubbard model 0,113. 
In the Hamiltonian, we have already used the LDA, i.e., 
u v (x) = u v (n(x)) and K v (x) = K v (n(x)). The canon- 
ical momenta U u are conjugate to the phases <j) Uj i.e., 
[(j) v (x) , U p (V)] = iSpyS (x — x'). Physically, the gradi- 
ents of the phases d x (j) u (x) are proportional to the den- 
sity (or spin density) fluctuations, Sn u (x) = —d x (\) v (x), 
while the momenta U u (x) are proportional to the density 
(or spin density) currents, j u (x) = u u (x)K u (x)U u (x). 

From this linearized Hamiltonian, we derive the hydro- 
dynamic equation of motion: 



d 2 5n v 
Or 2 



_d_ 

dx 



u u (x)K u (x)- 



d ( u v [x) 



dx \K u (x) 



5n v 



(5) 



Hereafter we focus on the readily observable density 
modes. We consider the n-th eigenmode with Sn p (x) ~ 
Sn p (x) exp {iuo n r) and substitute it into Eq. The 
resulting equation for Sn p (x) is solved using boun dary 
conditions which require the current j p (x) to vanish |ll| . 
The details, together with the interesting question of 
spin-charge separation, will be presented elsewhere (T^ . 

Fig. 2 shows the square of frequencies of the breathing 
mode and of the dipole mode at fixed values of v and of k. 
For a constant value of ft, the remarkable feature shown 
in Figs. 2b and 2d is the appearance of a sharp dip at a 
particular value of the characteristic filling factor, which 
has been identified as the transition point from a metal 
to an insulator. On the left side of the transition point, 
the cloud is a pure metal and the density fluctuates at 
all sites. The decrease of the frequency with increasing v 
can be understood from the fact that in the homogeneous 
case the cloud has a sound-like spectrum, and the sound 
velocity u p decreases as u p oc 1 — n on approaching the 
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Figure 2: Square of frequencies of the breathing mode (Figs. 
2a and 2b) and of the dipole mode (Figs. 2c and 2d) at fixed 
values of v and of n . The circles in (a) are the predictions of 
the sum-rule approach, while the circles in (b) and (d) are the 
results for the non-interacting Fermi gas, obtained by solving 
the semi-classical Boltzmann kinetic equations. 



transition point at n = 1, due to the opening of the 
charge gap. For larger filling, an insulating domain forms 
at the center of the trap, in which the LL description 
breaks down. In response to a trap displacement, the 
center of mass will develop a quasi-static displacement 
in this case. Our theory is able to describe the local 
oscillations of the metallic wings around this new origin. 

We have verified this expectation by using the semi- 
classical Boltzmann kinetic theory for a non-interacting 
cloud As v increases, the metallic wings shrink in 

size and become less compressible. Therefore the mode 
frequency rises after crossing the transition point. This 
is the central result of this Letter. The Boltzmann the- 
ory also lets us investigate nonlinear damping effects, 
which we find increase rapidly with trap displacement 
(T^ . These are absent in the LL theory, which only treats 
small displacements. 

The application of the LL theory for low-energy dy- 
namics of the ID Hubbard model in the unconfined case 
is already well accepted in condensed matter physics fTc| . 
In the presence of a trap, we would expect that the LDA 
is still applicable for describing the dynamics of the cloud. 
To this aim, we have examined the validity of the LL 
theory within LDA for collective oscillations by using al- 
ternative methods in two limiting cases, (i) In Fig. 2a 
the open circles are the prediction of m^/mi or mi/m_i 
sum-rules at v — > [l^ 0|. The excellent agreement 
shows that the LL theory is correct at low density for ar- 
bitrary coupling strengths, (ii) On the other hand, for a 
non-interacting gas at a finite filling factor, we calculate 
the mode frequency by solving the semi-classical Boltz- 
mann kinetic equation pjj. As shown by open circles 
in Figs. 2b and 2d, this leads again to the same result 
obtained by the LL theory. Combining these two obser- 
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Figure 3: (color online) A proposed experimental scheme. 
The square of breathing mode frequency has been plotted 
as a function of the depth of the lattice potential s for a gas 
of potassium atoms. The inset shows the path in the phase 
diagram with increasing s. 

vations, we expect that the LL theory within the LDA is 
valid for finite filling factor at arbitrary couplings, under 
the required conditions. 

Experimental proposal — In recent experiments, an 
array of quantum gases of bosonic 87 Rb atoms in ID op- 
tical lattices has been created [15]. The gas is first sub- 
jected to a combined potential of a harmonic trap and 
deep transverse standing waves. An additional optical 
lattice with adjustable depth is then switched on in the 
axial direction to form the final configuration. An in- 
teracting Fermi gas of 40 K atoms in 3D lattices has also 
been demonstrated [ifl]. We expect that an array of in- 
teracting gases of 40 K atoms will be realized soon. Using 
the realistic parameters in Ref. [15J (with 87 Rb being 
replaced by 40 K ), we have calculated th, u = U/th, and 
u in a Gaussian approximation [13] . We have taken an 
inter-species s- wave scattering length asD = +26a# for 
40 K atoms, where a# = 0.529 A is the Bohr radius. The 
value of CL3D is tunable by a Feshbach resonance. The 
typical number of atoms in each tube is TV = 15, which 
is about the minimum required for the validity of LDA. 

We show the square of frequencies of the breathing 



mode as a function of the depth s of the lattice potential 
in Fig. 3. Finite temperatures and/or nonlinear effects 
due to finite amplitude oscillations, will lead to the de- 
cay of collective modes, and simultaneously weaken the 
sharpness of the frequency dip. In addition, the exper- 
iments currently average over many tubes of ID inter- 
acting gases, where the number of atoms has a parabolic 
distribution. This will further reduce the dip signal. Nev- 
ertheless, the qualitative picture of the appearance of dip 
structure in the vicinity of MMIT transition should per- 
sist, and will be a signature of the entrance into the Mott- 
insulator phase. 

Conclusions. — In summary, we present a phase dia- 
gram for ID confined fermions in an optical lattice. We 
show that the behavior of the frequencies of collective 
density modes is distinct in each phase. Therefore, the 
phase diagram can be detected unambiguously by mea- 
suring density oscillations. This provides a useful tool for 
locating the quantum phase transition from the metallic 
phase to the Mott-insulator phase. We expect that a 
similar sensitivity of the mode frequencies with respect 
to the phase will arise in higher dimensions. Further- 
more, this signature could be applied as well to a gas of 
ID bosons in an optical lattice, since the behavior of its 
low-energy excitations in the superfluid phase falls in the 
same universality class as the LL theory. 

Our investigation is based on the exact Bethe ansatz 
solution and the LL theory of ID Hubbard model, with 
the effect of harmonic traps being incorporated in the 
LDA. We justified the use of LDA in several limiting 
cases, both for static and dynamical properties. To fur- 
ther examine the validity of LDA for dynamics, more 
complicated numerical simulations are necessary. Pos- 
sible candidates are fermionic phase-space [18] or time- 
dependent density-matrix renormalization group method 
(lflj . This is beyond the scope of the present Letter. 

We are indebted to M. Rigol and F. Werner for use- 
ful discussions. X.-J. L and P. D. D gratefully acknowl- 
edge support by the Australian Research Council. H. 
H was partially supported by INFM under the PRA- 
Photonmatter Programme. 
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